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We have computed phase diagrams for rotating spin-1 Bose-Einstein condensates with long-range 
magnetic dipole-dipole interactions. Spin textures including vortex sheets, staggered half-quantum- 
and skyrmion vortex lattices and higher order topological defects have been found. These systems 
exhibit both superfluidity and magnetic crystalline ordering and they could be realized experimen- 
tally by imparting angular momentum in the condensate. 

PACS numbers: 03.75.Lm, 67.85.De 
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A bar magnet is an omnipresent dipolar object in clas- 
sical physics. By dislodging one magnet in the proximity 
of another the striking dependence of the interaction be- 
tween dipoles on their separation and relative orientation 
becomes apparent. In addition to magnetism, dipolar ef- 
fects play role in a multitude of physical systems and pro- 
cesses including wireless telecommunications, spintron- 
ics, chemical bonding between atoms and molecules, and 
liquid crystals pQ. 

In the recent years, several candidates exhibiting dipo- 
lar effects in ultra-cold quantum gas experiments have 
emerged [2J. Polar molecules with permanent electric 
dipole moments [3] may provide a system amenable to 
a semiclassical treatment 01 [5] . Chromium condensates 
have been observed to show dipolar effects due to their 
intrinsically large magnetic moments [6l [7] with dyspro- 
sium [5] and erbium promising even stronger dipolar ef- 
fects and elevated potential for quantum simulators and 
quantum information processing applications [2] ■ It may 
also be possible to achieve co-existence of both superfluid 
and magnetic long-range order in these systems and to re- 
alize elusive supersolid states — hints of which has already 
been observed in experiments with spinor 87 Rb [9J[TU]. In 
addition, dipolar cold atomic clouds arc predicted to fea- 
ture a range of quantum phenomena previously unseen 
in Bose-Einstein condensed gases [TTUI4] . 

Here we calculate phase diagrams for rotating spin-1 
Bose gas in the presence of long-range magnetic dipole- 
dipole interactions without constraining the spin degree 
of freedom. On changing the relative strength of the long- 
range dipole-dipole interaction with respect to the local 
contact interaction, structural lattice transitions occur in 
the system. Due to the dipolar interactions a variety of 
textures including vortex sheets, staggered skyrmion and 
half-quantum vortex lattices and higher order topological 
defects emerge in these systems. 

Model — We consider a zero-temperature spin-1 Bose- 
Einstein condensate and model it using a mean-field the- 
ory which includes magnetic dipole-dipole interactions. 
Such system is described by a three-component spinor 



order parameter $ = ["0i(r), tpo( r ), ' ( /'-i( r )] T [T5I [TB], 
The ground-state wavefunction is a solution to a dipo- 
lar spinor Gross-Pitaevskii equation 
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where the operator £(r) = — ft 2 V 2 /2m + K x t(r) — QL Z — 
/i + gn(v) contains a single particle part and a contact 
interaction term, where n(r) = J2k iV'fclr)! 2 is the total 
particle density. Notwithstanding there does not exist, 
due to dipolar interparticle potential, a rotating refer- 
ence frame where the laws of equilibrium statistical me- 
chanics would apply |17) . we may assign Q to be a La- 
grange multiplier associated with the angular momen- 
tum of the system measured by the operator L z . The 
external potential V r ext (r) = m[uj^_(x 2 + y 2 ) + u 2 z 2 ]/2 is 
fixed by the transverse uj± and axial cu z harmonic fre- 
quencies. The normalization condition of the wavefunc- 
tion, J v n(r)dr = N, relates the number of particles N 
in the system to the chemical potential /i. The spin- 
exchange interaction yields the terms f z (r) = g s (F z ) and 
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i(F y ))/y/2. The terms d z (r) = g d I z 



&ndd xy (r) 



iI y )/V2, where I k (r) = fdr'({F k )- 
, .. .... ,...,., i ■ -'| 3 ) are due to the dipole-dipole in- 
teraction. Here the vector e = (r — r')/|r — r'| and Fk 
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r — r')/|r — r' 
denote the components of the spin-1 operator F and (•) 
are their expectation values in the ground state. The cou- 
pling constants g = Airh 2 (an + 2a2)/3m, g s — 4irh 2 (a2 — 
ao)/3rn and ga — MoMlfff / '^ determine, respectively, 
the strength of s-wave, spin-exhange and dipole-dipole 
interactions between particles of mass m in terms of the 
bare scattering lengths a^ and a-i and the permeability 
of vacuum fiQ, Bohr magneton /zb and the Lande factor 

9F- 

We use hu± and a± = y^h/mui^ for the units of en- 
ergy and length, respectively. Our computational sys- 
tem is specified by the aspect ratio u z = 10 uj± and 
dimensionless parameters gN/fko±_a\ = 1000, g' s = g s j g 



0.6 
0.5 


(a)\ 


V 


(i,D 


0.4 
fio.3 


Illb \ 
(1,1) \ 


(1,0) 


■"ilia 


0.2 






(1,1) ' 


0.1 



1C 


1 (0,0) 


II (0,1) 


1C 


V 


10"' 



0.5 


(b) 




V 

(1,1) 


0.4 


IVa 


IVb 


(1,0) ^ 


fio.3 


QAJA) 


(«,%) 


^^Illa 


0.2 






(1,1) ' 


0.1 








I (0,0) 


II (0,1) 


1C 


- 3 


10" 


J 10"' 



9d 



FIG. 1. Schematic phase diagrams in the g' d — Q! parameter 
space for (a) ferromagnetic g' s = —0.01 and (b) antiferromag- 
netic g' s — 0.03 spin-exchange interaction. The numbers de- 
note the type of topological defects in each phase as described 
in the text. Phase boundaries are estimated by straight lines. 



and g' d = gd/g- As in [T8], we choose g' s = 0.03 for an- 
tiferromagnetic and g' s = —0.01 for ferromagnetic spin- 
exchange interaction corresponding to the values for 23 Na 
and Rb, respectively. The dimcnsionless frequency 
Q,' = Q/loj_ and g' d are the adjustable parameters of our 
numerical experiments. We assume that external mag- 
netic fields are absent and therefore exclude linear and 
quadratic Zeeman terms from Eq.(nj). 

Elementary vortex types — To facilitate further discus- 
sion we categorize certain topological defects relevant to 
our system. Invariance of the spinor order parameter \1/ 
under a suitable combination of spin rotation and gauge 
transformation W' = e j6 'e~ l ^ n ' F, I', where n is a unit vec- 
tor along the spin-rotation axis, allows us to classify the 
line defects in our system by associating with them a pair 
of winding numbers (£, s). Here £ — \9\/2tt measures the 
strength of mass circulation about the defect line and 
s — \c/)\/2ir is the spin winding number. With this nota- 
tion, a vortex which carries one unit of orbital angular 
momentum per particle and involves a 2tt spin rotation 
is denoted by (1,1). Such vortex has a ferromagnetic 
core and depending on the underlying model and imposed 
boundary conditions is known as a Mermin-Ho- |19j or 
Anderson- Toulouse vortex, [20] or a skyrmion [21] . These 
have recently been created in transition metal compound 
MnSi [53] and in BECs as dynamical states [53] ■ 

A half-quantum vortex or an "Alice string" is charac- 
terized by the doublet (5, §) and is predicted to exist e.g. 
in superconductors, superfiuid helium systems and spinor 
BECs pMH2l)] . The half-quantum vortex has a ferromag- 
netic core structure and vanishing spin density elsewhere. 
Half-quantum vortices have recently been created in po- 
lariton condensates [27] , 

In addition to the two "mixed" spin-mass textures, 
pure spin vortices (0,1) and pure mass vortices (1,0) 
emerge in our system. The former has a core devoid 
of spins and is therefore said to possess an unmagne- 
tized polar core, whereas the spin currents vanish around 
the latter. All of the above described topological objects 
are"coreless vortices" in the sense that the total particle 
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FIG. 2. (Color online) Skyrmion textures for g' B — —0.01, 
f2' = 0.6 and (a) g' d = 0.001 and (b) g' d = 0.01 corresponding 
to phases Illta and V, respectively. The main frames show 
the spin texture, arrows indicating the orientation of local 
spin and contour lines are plotted for the ^-component of 
spin density. Insets to the right of the main frames show the 
particle density n(r) and spin density \/^2 k {Fk} 2 - All square 
frames have the same (7a± x 7a±) spatial extent. 



density is finite at the vortex cores [2"5] . 

Ferromagnetic spin- exchange interaction — Figure [I] 
(a) shows a schematic phase diagram in the g' d — O' 
parameter space for a ferromagnetic g' s = —0.01 spin- 
exchange interaction. For O' = two phases exist de- 
pending on the relative strength of the dipolar interac- 
tions [U [12] [18] . In I the ground states are free of topo- 
logical defects whereas in II spontaneous spin currents 
emerge in the system. The ground states in II are occu- 
pied by a single (0,1) vortex [TB]. Eventually, for g' d s=s 1/4 
the system becomes unstable due to the effective attrac- 
tion caused by the dipole-dipole term [2] [5] [TH] [29] . The 
critical value of g' d is fairly insensitive to the values of $7' 
and g' s for the range of parameters studied here. In ex- 
periments, the value of g' d has been successfully tuned for 
spin-polarized chromium condensates [7] . In the absence 
of magnetic fields this could be achieved using optical 
Feshbach resonances [30] . 

For g' d = on increasing fi' above a critical value, mass 
currents are nucleated in the system [3T]. For ferromag- 
netic spin-exchange term g' s < the ground states in Ilia 
contain a single (1, 1) vortex. The phase Ilia precedes the 
dipolar collapse whenever Q' exceeds the critical value. 
For small values of g' d on increasing $7' further staggered 
skyrmion lattices composed of (1,1) vortices form as il- 
lustrated in Fig. [2la) which displays a typical phase Illb 
texture. The arrows in the main frame show the local 
spin orientation in the z — plane. The contour lines 
are plotted for the z-component of the local spin den- 
sity. The inset in the top right corner shows the particle 
density. Due to the weak dipolar coupling and since the 
skyrmion cores are filled with spins which are oriented 
parallel and antiparallcl to the z-axis, the spin density 
shown in the bottom right corner is rather smooth. The 
spin-texture in Fig. Wt&) is seen to be composed of two 
interleaved skyrmion lattices of opposite vortex core po- 



larisations and spin windings. Generally, individual (1,1) 
vortices do not possess cylindrical symmetry in these sys- 
tems. On further increasing the value of g' d keeping Of 
fixed, type (1,0) vortices enter the system [32]. A vari- 
ety of textures composed of a mixture of (1, 1) and (1, 0) 
vortex molecules are observed in phase V, an example of 
which is illustrated in Fig. [21(b). The (1,0) vortices leave 
a clear signature in the spin density due to their spinless 
core structure. 

Antiferromagnetic spin-exchange interaction — For an- 
tiferromagnetic spin-exchange term g' s > the phase di- 
agram shown in FigflTa) reduces to that of the ferromag- 
netic one shown in FigpTb) both for small fi' and for large 
g' d where dipolar effects dominate over the spin-spin in- 
teraction. In phase IVa the constituent defects are of type 
(h,h) and their spins lie predominantly in the x-y plane. 
The spins in the cores of the nearest neighbor defects are 
pointing in the opposite directions forming a staggered 
antiferromagnetic structure as shown in Fig(3[a). An im- 
age of the spin density strikingly reveals the underlying 
magnetic crystal. Magnetization density in these rotating 
states is analogous to the fractional 1/3— vortex lattice 
states found in the cyclic state of F — 2 spinor conden- 
sates [33] . On crossing from IVa to IVb the square lattice 
of (5,5) vortices is first squeezed into a rhombic shape 
due to the head-to-tail attraction of the dipoles. Next 
the staggered magnetic crystal phase IVa melts and the 
vortices rearrange to form a vortex sheet in phase IVb. In 
a vortex sheet the spin axes of the (5,5) vortices tend to 
align with the orbital angular momentum. Figure a\ (b) 
illustrates a typical texture emerging in phase IVb. Sim- 
ilar vortex sheets were first observed in superfluid 3 He-A 
[31]. For even larger values of g' d , in phase V, the {\,\) 
vortices are converted to (1,1) and (1,0) defects where 
the former (latter) can be thought to be constructed by 
merging two (|, \) vortices with same (opposite) spin 
windings. In the same spirit, the (0, 1) polar core vortices 
cannot be built from (| , |) vortices since the external ro- 
tation breaks the time reversal symmetry. 

Higher order topological dejects — Rapid rotation of 
the system leads to a centrifugal explosion of the cloud 
whereas strong enough dipole-dipole interaction causes 
the condensate to implode. Therefore the simultaneous 
limit Q! — > 1 and g' d — > 1/4 deserves particular atten- 
tion. Approaching this limit, Fig[4] shows textures under 
rapid rotation in strongly dipolar systems. In Figffl (a), 
we interpret each loop of the 'figure of eight' as deformed 
type (4, 1) defect and in FigJ4] (b) a single (5, 1) defect 
is identified. In addition, we have found a wide variety 
of both axisymmetric and non-axisymmetric high wind- 
ing number defects. It is thus conceivable that a whole 
hierarchy of higher order topological defects emerge in 
this dipolar system. In the absence of dipolar forces, 
higher order topological objects are suppressed due to 
the quadratic increase of the kinetic energy as a func- 
tion of the multiplicity of the topological defect. Here 
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FIG. 3. (Color online) Magnetic crystal and vortex sheet 
textures respectively for (a) g' d — 0.001 and (b) g' d — 0.01 
corresponding to phases IVa and IVb. Contour lines in (a) and 
(b) are plotted for the x- and z-components of spin density, 
respectively. For both (a) and (b) g' s — 0.03 and il' — 0.6. 
All square frames have the same (7a± x 7a±) spatial extent. 
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FIG. 4. (Color online) Higher order topological defects cal- 
culated for (a) g' d = 0.23 and (b) g' d = 0.24. For both 
states g' s — —0.01, Q! = 0.95, contour lines are plotted for 
the z-component of spin density and all square frames have 
(5a± x 5a±) spatial extent. 



the addition of dipolar interactions enables the existence 
of such textures irrespective of the sign of the spin-spin 
coupling. Furthermore, due to the intrinsic effective con- 
finement provided by the dipolar interactions which, to- 
gether with the trapping potential, may compensate for 
the centrifugal expansion, the experimental realization of 
the lowest Landau level physics may become reachable in 
these rapidly rotating dipolar systems. 

In conclusion, we have calculated phase diagrams for 
rotating dipolar F — 1 spinor Bose-Einstein condensates. 
For weak dipolar coupling and for ferromagnetic (antifer- 
romagnetic) spin-exchange interaction the spin textures 
composed of skyrmions (half-quantum vortices) arrange 
in an interlaced double square lattice structures where 
the magnetized vortex cores exhibit staggered antifer- 
romagnetic ordering. The orbital angular momentum 
loaded in the system causes quantized mass currents to 
nucleate in the condensate due to its inherent superfluid 
property. As a consequence magnetic crystalline order- 
ing is induced in the system by the spin vortices gener- 
ated via the spin-gauge symmetry. These steadily rotat- 



ing systems are superfluid and form equilibrium magnetic 
crystals. They also exhibit spatial periodic modulation in 
the particle density. For antiferromagnetic spin-exchange 
coupling on increasing the dipolar interaction strength 
the magnetic crystals melt and form vortex sheet tex- 
tures. Various spin structures emerge already for small 
values of relative dipolar interaction strength and might 
therefore be observable experimentally. Our results, al- 
beit calculated for rotating systems, may also shed light 
to the interpretation of the short-range magnetic struc- 
tures observed in [5], 

Lattice defects such as interstitials and vacancies fre- 
quently present in the non-ground states leave an ex- 
tremely clear signature in the particle density profiles 
and could be detected using Stern-Gerlach type ex- 
periments in combination with standard imaging tech- 
niques. Therefore these systems might also provide an 
excellent platform to investigate proliferation of defects 
in the finite-temperature Berezinskii-Kosterlitz-Thouless 
superfluid transition [55H4T] . 

For strong enough dipole-dipole coupling the system 
collapses due to the induced effective attraction whereas 
rotation exceeding the harmonic trapping frequency leads 
to centrifugal explosion. This suggests the possibility 
of entering exotic quantum-Hall like states for which 
the dipolar attraction and centrifugal effect would be 
counter-balanced. Approaching such limit we have found 
a wide variety of higher order topological defects entering 
the system. 

This work was supported by the Japan Society for the 
Promotion of Science (JSPS). 



[1 
[2 
[3 

K 

[5 

[7: 

[8 

[9 

[10 



P. G. de Gennes and J. Prost, The Physics of Liquid 

Crystals, (Oxford University Press, Oxford 1993) 

T. Lahaye, C Menotti, L. Santos, M. Lewenstein, and T. 

Pfau , Rep. Prog. Phys. 72, 126401 (2009). 

K. Ni, S. Ospelkaus, M. de Miranda, A. Pe'er, B. Neyen- 

huis, J. Zirbel, S. Kotochigova, P. Julienne, D. Jin, and 

J. Ye, Science 322, 231 (2008). 

M. Takahashi, S. Ghosh, T. Mizushima, and K. Machida, 

Phys. Rev. Lett. 98, 260403 (2007). 

J. A. M. Huhtam&ki, M. Takahashi, T. P. Simula, T. 

Mizushima, and K. Machida, Phys. Rev. A 81, 063623 

(2010). 

J. Stumer, A. Griesmaier, T. Koch, M. Fattori, T. Pfau, 

S. Giovanazzi, P. Pedri, and L. Santos, Phys. Rev. Lett. 

95, 150406 (2005). 

T. Lahaye, J. Metz, B. Frohlich, T. Koch, M. Meister, 

A. Griesmaier, T. Pfau, H. Saito, Y. Kawaguchi, and M. 

Ueda, Phys. Rev. Lett. 101, 080401 (2008). 

M. Lu, S. H. Youn, and B. L. Lev, Phys. Rev. Lett. 104, 

063001 (2010). 

M. Vengalattore, S. R. Leslie, J. Guzman, and D. M. 

Stamper-Kurn, Phys. Rev. Lett. 100, 170403 (2008). 

M. Vengalattore, J. Guzman, S. Leslie, F. Serwane, and 

D. M. Stamper-Kurn, arXiv:0901.3800 



[11 
[12 

[13 

[14; 

[15 

[16 
[17 

[18 
[19 

[20 

[21 

[22 

[23 
[24 
[25 
[26 
[27 

[28 
[29 

[30; 

[31 
[32 
[33 
[34 

[35 

[36 
[37 
[38 
[39 
[40 
[41 



L. Santos, G. V. Shlyapnikov, and M. Lewenstein, Phys. 

Rev. Lett. 90, 250403 (2003). 

Y. Kawaguchi, H. Saito and M. Ueda, Phys. Rev. Lett 

97, 130404 (2006). 

Y. Kawaguchi, H. Saito and M. Ueda, Phys. Rev. Lett 

96, 080405 (2006). 

L. Santos, T. Pfau, Phys. Rev. Lett. 96, 190404 (2006). 

T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 

(1998). 

T. L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

A. J. Leggett, Quantum Liquids (Oxford University 

Press, Oxford 2006). 

S. Yi and H. Pu, Phys. Rev. Lett. 97, 020401 (2006). 

N.D. Mermin and T.-L. Ho, Phys. Rev. Lett. 36, 594 

(1976). 

P. W. Anderson and G. Toulouse, Phys. Rev. Lett. 38, 

508 (1977). 

T. H. R. Skyrme, Proc. R. Soc. Lond. A 260, 127 (1961); 

Nucl. Phys. 31, 556 (1962). 

A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, R. Ritz, 

P.G. Niklowitz, and P. Boni, Phys. Rev. Lett. 102, 

186602 (2009). 

L. S. Leslie, A. Hansen, K. C. Wright, B. M. Deutsch, 

and N. P. Bigelow, Phys. Rev. Lett. 103, 250401 (2009). 

G. E. Volovik, The Universe in a Helium Droplet 

(Clarendon Press, Oxford, 2003). 

T. Isoshima and K. Machida, Phys. Rev. A 66, 023602 

(2002). 

J. Ruostekoski and J.R. Anglin, Phys. Rev. Lett. 91, 

190402 (2003). 

K. G. Lagoudakis, T. Ostatnicky, A. V. Kavokin, Y. G. 

Rubo, R. Andre B. Deveaud-Pledran, Science 326, 974 

(2009). 

The (1,0) vortices in our system are nonaxisymmetric 

"split-core" states having finite core-density. 

N. G. Parker, C. Ticknor, A. M. Martin, and D. H. J. 

O'Dell, Phys. Rev. A 79, 013617 (2009). 

M. Theis, G. Thalhammer, K. Winkler, M. Hellwig, G. 

Ruff, R. Grimm, and J. H. Denschlag, Phys. Rev. Lett. 

93, 123001 (2004). 

T. Mizushima, N. Kobayashi and K. Machida, Phys. Rev. 

A 70, 043613 (2004). 

T. Mizushima, K. Machida and T. Kita, Phys. Rev. Lett. 

89, 030401 (2002). 

J. A. M. Huhtamaki, T. P. Simula, M. Kobayashi, and 

K. Machida, Phys. Rev. A 80, 051601(R) (2009). 

\j. Parts, E. V. Thuneberg, G. E. Volovik, J. H. Koivu- 

niemi, V. M. H. Ruutu, M. Heinila, J. M. Karimaki, and 

M. Krusius, Phys. Rev. Lett. 72, 3839 (1994). 

Z. Hadzibabic, Z. Kriiger, M. Cheneau, B. Battelier, and 

J. Dalibard, Nature 441, 1118 (2006). 

P. Clade, C. Ryu, A. Ramanathan, K. Helmerson, and 

W. D. Phillips, Phys. Rev. Lett. 102, 170401 (2009). 

T. P. Simula and P. B. Blakie, Phys. Rev. Lett. 96, 

020404 (2006). 

T. P. Simula, M. J. Davis, and P. B. Blakie, Phys. Rev. 

A 77, 023618 (2008). 

R. N. Bisset, M. J. Davis, T. P. Simula, and P. B. Blakie, 

Phys. Rev. A 79, 033626 (2009). 

M. Holzmann and W. Krauth, Phys. Rev. Lett. 100, 

190402 (2008). 

V. Pietila, T. P. Simula, and M. Mottonen, Phys. Rev. 

A 81, 033616 (2010). 



